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In this letter we suggest that the obvious discrepancy between theoretical prediction on the non - DD decays 
of 1^(3770) and data is to be alleviated by taking final state interaction (FSl) into account. By assuming that 
- - - 1^(3770) overwhelmingly dissociates into DD, then the final state interaction induces a secondary process, we 

calculate the branching ratios of i/'(3770) — > DD — > J/if/r], pn, air], K'K. Our results show that the branching 
O , ratio of i/r(3770) -> non - DD can reach up to S™^^ = (0.2 ~ 1.1)% while typical parameters / = 0.4 GeV"^ 

' and a = 0.8 ~ 1.3 are adopted. This indicates that the FSl is obviously non-negligible. 

in ! PACS numbers: 13.30.Eg, 13.75.Lb, 14.40.Lb 

Ph. 

Obviously, physics about Charmonium is still not a closed book yet, instead, this field is full with challenges and opportunities 
■ ill]- Due to developments and improvements of facility and technique, the detection precision is greatly enhanced in the past 
decade, as a direct consequence new puzzles have continually emerged. Indeed, some old puzzles have been understood now, 
I— I but a number of them remain unsolved yet. Theorists are endeavoring to look for solutions. The general strategy is that first, one 
can fumble solutions in the framework of present theory i.e. QCD and see if we miss something in our calculations, then if all 
possibilities are exhausted one needs to invoke new physics beyond the standard model. In this work, we follow the first strategy 
to explain excessive non-DZ) component of inclusive i/'(3770) decay, which is conducted in a series of experiments by the BES 
1^ collaboration 13. H 0. S 0] in the past three years. 

As a well measured charmonium state, i^(3770) generally is considered as a mixture of 2^^5'i and l^Di states HE]- Since 
3770 MeV is a bit above the threshold of DD production, such a bound state may dissolve into open charms which eventually 
hadronize into DD. Therefore, before observing sizable non-DD decay rates, 11/(3770) was supposed to dominantly decay into 
> DD, including and D^D . There could be some possible non-DD modes fl [M [U E [H [P], especially the hidden 

O charm decay modes, such as J/if/jm and J/i/zj] with S[(A(3770) ^ J/if/n+Ti-]^ (1.93 + 0.28) x 10-^ S[iA(3770) ^ Jlipn^n^] = 
00 (8.0 ± 3.0) X 10^"* and S[i/r(3770) J/il/rj] = (9 + 4) x lO"'* respectivelyjpj], and El radiative decays yxcj with decay widths 
1 72 ± 30 keV, 70 + 1 7 ke V and < 2 1 keV for 7 = 0, 1,2 respectively iH M. The sum of all the branching ratios of these hidden 
charm decay modes is less than 2%, so all these measurements support the allegation that i/'(3770) overwhelmingly decays into 
(N DD. 

However, the BES collaboration investigated the inclusive decays of i/r(3770) and found that the branching ratio of i^(3770) — > 
^ DD is about (85+5)% HH. This is later verified by the measurements of non-DD inclusive processes with the branching fraction 
S[iA(3770) non - DD] = (13.4 ± 5.0 + 3.6)% |5] and S[tfr(3770) non - DD] = (15.1 ± 5.6 ± 1.8)% |6] respectively by 
adopting two different methods. The CLEO measurements indicate cr(e+e" — > i/r(3770) — > hadrons) - (6.38 + 0.08^q3q) nb |lg] 
andcr(e+e- ^ (/r(3770) ^ DD) = (6.57 ±0.04 ±0.0 l)nb id, which together make a S(iA(3770) ^ DD) = (103.0 ± 1.4+5 ')%. 
5J] Notice that the error on the high side is about 6.8% [37], by this error tolerance, there could be a large (10 ~ 15)% fraction of 
i/'(3770) non - DD decays. The CLEO and BES results are inconsistent at > 2cr level, and we would employ both of them as 
inputs to our numerical computations and an obvious difference is explicitly noticed. 

Kuang and Yan 1 18] calculated the i^(3770) — > J/ifrnn using the QCD multi-expansion, which properly deals with the emission 
of light hadrons during heavy quarkonia transitions (for a review see an enlightening paper 1 19]). Their prediction is consistent 
with the exclusive measurement on hidden charm decays of i^(3770). It is generally concuiTed that, the measurements on 
the well measured channels J/ifmn, J/iffrj and yXcj are consistent with present theoretical predictions. Thus to understand the 
experimental results, one should find where (^(3770) goes besides J/tj/nn, J/iffrj and JXcj- Recently He, Fan and Chao ll20ll 
introduced the color-octet mechanism and calculated the i^(3770) — > light hadrons in the framework of NRQCD by considering 
next to leading order contribution. The calculation result shows that r[i/'(3770) — > light hadrons] is 4-(i7^l^l keV. If combing 
radiative decay contribution with that of (^(3770) — > Light Hadrons, the branching ratio of the non-DD of i/r(3770) is about 5% 
ll20ll . which is still three times smaller than 15% non-DD branching ratio measured by the experiment. 
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Instead, Voloshin suggested, ^{3110) is not a pure cc state. There exists a sizable four quark component (mm + dd)cc and the 
fraction is about (9(10%) in i//{3nO), which results in a measur able rate of iA(3770) ^ n^J/ip, t]J/i// |21]. 

Generally, one can categorize the strong decay modes of i/r(3770) into three types: open charm decay (DD), hidden charm 
decay J/iJ/X {X - light mesons) and the decay into light hadrons (L-H decay). One can be more confident that the rates of hidden 
charm decays are properly evaluated in terms of the QCD multi-expansion, and the L-H decay occurs via three-gluon emission 
mechanism cc — > 3g. 

There is an alternative explanation to the puzzle. Twenty years ago, Lipkin proposed that the non-DD strong decays of (/'(3770) 
realize via DD intermediate states, and further suggested that 1/^(3770) does not 100% decay into DD |22]. Later Achasov and 
Kozhevnikov calculated the non-DD channels of (/'(3770) only considering the contribution from the imaginary part of the decay 
amplitude |23]. Namely such final state interactions which are involved in the hadronic loop effects, do contribute to both the 
hidden charm and L-H decays. The essential point of the loop effect is attributed to the coupled channel effects. A quark-level 
process is explicitly illustrated in the left diagram of Fig. [T] Such a mechanism should exist in all hidden charm and L-H decays 
of charmonia ll24l 125 1 . As shown in Fig. [T] iff{3770) — > J/i//X and i/'(3770) — > Light Hadrons processes do not suffer from the 
Okubo-Zweig-Iizuka (OZI) suppression. Since 1^/(3770) DD takes place near the energy threshold, one can expect that the 
FSI may be significant. 




FIG. 1: Quark-level descriptions of the hadronic loop mechanism for the hidden charm decay (diagram (a)) and L-H decay (diagram (b)) of 
(^(3770). 

In this letter, we focus on two-body hidden charm decay modes (J/il/ri) and two-body L-H decay modes (pn, ojj] and K*K) 
which obviously are the main ones. Here K*K denotes K*K + K*K. 

In order to calculate the hadronic loop effect in strong decays of i/'(3770), we consider the diagrams depicted by Fig. |2] 
which are an alternative description in the hadron-level language. 1/^(3770) first dissolves into two charmed mesons, then by 
exchanging D* in t-channel, they turn into two on-shell real hadrons ^ and S. Since the dissociation does not suffer from the 
OZI suppression, one can expect it to be dominant. 




FIG. 2: The hadron-level diagrams depicting hadronic loop effect on i/r(3770) — » — » J/if/rj, p^n". Of course they can be simply replaced 
by other states D'*>''(Z)<*"') -> D<*'+(£)(*'-) with an isospin transformation and D**' — » D'*' with a charge conjugate transformation to constitute 
new but similar diagrams. By replacing relevant mesons, we can obtain the diagrams for a)r] and K'K channels. 



One can obtain the absorptive part of the decay amplitude of 1/^(3770) -^D + D^^ + S - J/i//?], pn, ojt] and K*K) 

AmeiM]) = J d^M* [(^(3770) ^ DD\ 

xM [dD ^ yis] ■ ^ (J) 

where |p| - [/l(M^, w^, OT^)j ^ /(2M^) is the three-momentum of the intermediate charmed mesons in the center of mass frame 
of i/'(3770). A(a, b, c) = + + c^ - lab - lac — Ibc is the Kallen function. The form factor T [w^,,,, <7^j is the key point for 
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the evaluation of the amphtude. One can use the monopole form factor (FF) 

rHA-^P^ (2) 

which compensates the off-shell effect of exchanged meson and describes the structure effect of the interaction vertex. As a free 
parameter, A can be parameterized as A(m,) = w,- + qKqcd [26] . m, denotes the mass of exchanged meson, Aqcd = 220 MeV. 
The range of dimensionless phenomenological parameter a is around 0.8 < a < 2.2 ll26ll . As a matter of fact, there are other 
possible forms for T [m^,,,, <7^j, such as the exponential one etc. in literature. Generally they are equivalent somehow, as long 
as their asymptotic behaviors are the same. 

Since the mass of i/'(3770) is close to the threshold of DD production, the dispersive part of the amplitude of iliiillQ) — » 
D + D — » y[ + S makes a large contribution to the decay width. By unitarity, one can obtain the dispersive part in terms of the 
dispersion relation. The total decay amplitude of i/'(3770) D + D ^ J[ + S which includes both absorptive and dispersive 
paits is expressed by ll27ll28ll29ll 



M[(/r(3770) -^DD^J{S\ 



_ 1 r 



—-—^dr+A:^s(M^, 



dr+Aj,siMi), (3) 



where r^^ - 4m^. After replacing in the amplitude in eq. ([TJ with r, we get the amplitude Aj3s(r). The energy dependent 

factor '??(r) is defined as 'R(r) - exp(-J |q(r)p), which not only reflects the |q(r)|-dependenceof the interaction between i/r(3770) 
and DD mesons, but also plays the role of ultraviolet cutoff. Meanwhile, '??(r) can be understood as the coupled channel effect 
summing up all the bubbles from the charmed meson loops ll30ll . Here |q(r)| denotes the three momentum of D meson in the 
rest frame of (/'(3770) with the mass x ^/r. The interaction length factor I is related to the radius of the interaction by 
I = R/6 1,31,1 . Pennington and Wilson indicated that I - 0.4 GeV"^ corresponding to R - 0.3 fm is favorable when studying 
the charmonium mass shift lIsTll . 

Based on the effective Lagrangian approach, one can formulate A1*[i/'(3770) — > DD] and M[DD — > ^I^SV The effective 
Lagrangians related to our calculation are constructed by considering the chiral and heavy quark symmetries ll32l[33l[34ll 

■Cd'dy = -2/».ov£^v./j id>'W% (if^IT"^ - D*ft&^ , 

where *P denotes charmonium states Jjijj and P and V are the octet pseudoscalar and nonet vector meson matrices, 

respectively. The values of coupling constants relevant to our calculation are g^oo - 4.70, gjj^-Dv - 4.25 GeV"', gjym - 
17.31 and/o-ov = 2.33 GeV-' determined in Refs. IMIlllll. 
For the process (/'(3770) — » D{k\) + D{k2) — > Jl^{k{) + 77(^4) by exchanging D*", one formulates its amplitude 
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cF - ml. 



The absorptive amplitude of (/r(3770) ^ D(ki) + D{k2) P(k3) + ¥(^4) reads as 



-8''' + ^]r'[ml.,q']. (4) 



X 



?M„. ; q^ - m^). 



Here the isospin factor from P and V matrices results in an extra factor Qjyb in the above amplitudes, which are 1 / V6, 1 /2, 1 
and 1/ VT2 for the amplitudes of J|^p^, p"7r°, K*^K- (K*-K+, K'^K'^, ^*°/:"), corj modes, respectively. If considering SU(3) 
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symmetry, the factor Qjis comes from both the isospin transformation D°^*^ £)(*)± and the charge conjugate transformation 
£)(•) ^ oi*)^ which results in (3y,s = 4, 2 for the ampHtudes of J/ifrrj (p^n^, corj) and K*^K- (K*-K^, K*'^K°, ^*"/:°) channels, 
respectively. 
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FIG. 3: The dependence of decay width of sum of i^(3770) — > i/i/'T/, pn, a>rj, K'K channels on a with several typical value I = 0.4 ~ 1.0 
GeV"^. In right-hand diagram, with I - 0.4 GeV"^, we also show the variation of partial decay widths of if/(3nO) — > J/if/rj, pn, ojt], K'K 
with a. 

In the left diagram of Fig. [3] we plot the dependence of numerical result on a with several typical value I = 0.4 ~ 1 .0 Ge V"^, 
which is the sum of decay widths of i^(3770) — > pn, corj, K*K induced by DD intermediate states. We set r[i/r(3370) — » DD 
pn] ^ 3r[i/r(3370) ^ DD ^ p'V'] and r[(A(3370) DD ^ K*K] ^ 4r[(A(3370) DD ^ K*+K-] which are determined by 
the SU(3) symmetry. The dependence of decay widths of each modes i^(3770) — » J/t/zT], pn, ajj], K*K induced by long-distant 
contribution on the parameter a within the range of 0.8 < a < 2.2 is shown in the right diagram of Fig. [3] where a typical value 
I = 0.4 GeV 2 is adopted. 

Our numerical results indicate that the decay widths of i^(3770) — > DD — > pn, K*K are about one order larger than that of 
i/'(3770) — > DD — > J/iffT], OJT] as we set the same values of parameters a and I for all the processes. The difference between the 
widths of tfr(311Q) — > DD — > pn, K*K and that of (^(3770) — » DD — > J/i/zt], ojt] is due to the phase space, factors @:^'B ™d Qy\'s- 
Whereas the amplitudes for pn, K'K are comparable, and they are the main non - DD decay channels obviously. 

The branching ratio of i/f(3770) non - DD including all J/iffJ], pn, corj, K*K modes with a fixed value I = 0.4 GeV"^ 
is shown in Fig. |4] within the range of a = 0.8 ~ 2.2. Furthermore, let us compare our result with the BES data 15i] and the 
result of iJ/{37TQ) Light Hadrons including the color-octet mechanism calculated up to next to leading order in the approach 
of NRQCD |20]. Fig. |4] shows that when FSI effects are taken into account, the NRQCD results plus FSI contribution can be 
very close to the BES data as long as a takes a value of 2.0~2.2. Since our results heavily depend on the parameter a, which 
is fully determined by the non-perturbative QCD effects and therefore cannot be determined based on a first principle, one can 
only phenomenologically fix it by fitting data. We also notice that the amplitudes in eqs. (Ill-Q are dependent on the values 
of coupling constant in every vertex, which results in that the decay width is proportional to the square product of all of the 
coupling constants. If the uncertainty is 20% for each coupling constant, the maximum of uncertainty the decay width is 4%. 



The BES data S[i/r(3770) pn] < 2.4 x 10"^ with corresponding width r[(/r(3770) pn] < 65 keV and the CLEO data 
S[if/(3110) pn] < 4.0 X lO""* with corresponding width r[^(3770) pn] < 109 keV (l2] help to further constrain the range 
of a to 0.8 < a < 1.1 and 0.8 < a < 1.3, respectively. The relevant values of the decay widths and the branching fraction are 
listed in TableU It is noted that when the FSI is taken into account and a is much restricted, the prediction of the branching ratio 
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FIG. 4: The comparison of our result (blue dash-dotted line) with the BES data (dashed line with shadow band) of excessive non-D5 component 
of the inclusive ifi(3770) decay |5] and the result of i/'(3770) — > Light Hadrons (dash-dotted line with shadowed band) by the color-octet 
mechanism calculated up to next to leading order within the framework of NRQCD j|20ll . Here the red line with green shadowed band is the 
total result including the NLO NRQCD effects and FSI contribution. The green shadowed band corresponds to the error tolerance, coming 
from the NRQCD estimate in Ref. i20ll . The orange and light blue shadowed bands are the suitable window for a, which is respectively 
determined by the BES data fs?] and CELO data 1 12] of i/^(3770) ^ pn. 



a 


Tpn (keV) 


r''^' „^ (keV) 

non-DD 


S"' (%) 

nnn-DD 


SINRQCD + Ours] (%) 


0.8 


20 


48 


0.2 


1.4 ~ 3.7(2.2) 


0.9 


32 


75 


0.3 


1.5 ~ 3.8(2.3) 


1.0 


47 


113 


0.4 


1.6 ~ 3.9(2.4) 


1.1 


66 


160 


0.6 


1.8-4.1(2.5) 


1.2 


94 


223 


0.8 


2.0 ~ 4.3(2.8) 


1.3 


127 


301 


1.1 


2.3-4.6(3.1) 



TABLE I: The typical values of decay width of pn channel r„„, the sum of decay widths F'^*' „, and the branching ratio S'^*' „^ of all channels 
discussed in this work. Here we fix / = 0.4 GeV"^. The branching fraction S[NRQCD + Ours], which is the sum of our result and the NRQCD 
result. The values in bracket are the central values. 



of (^(3770) non - DD caused by the FSI can reach up to S^^^^q = (0.2 ~ 1.1)% (taking CLEO data of i/r(3770) pn to 
constrain a). It indicates that even though FSI is significant, it cannot make a drastic change as long as a is restricted to be less 
than 2. 1 . Furthermore, the upper limit of the total contribution of the NRQCD and FSI is up to 4.6%. The branching ratios of El 
transition ip{311Q) yxcj (J=0,l,2) and ip{311Q) Jlipnn, Jjipr] are about_(1.5 ~ 1.8)% ll^iillill- If summing up all the 
above non-DD contributions, the branching ratio of the channels with non-DD final states can be as large as 6.4%, which is still 
smaller than the experimental value S[i/'(3770) — > non - DD] - (13.4 ± 5.0 + 3.6)% but near its lower bound 

As a short summary, let us emphasize a few points. First, even including contributions of color-octet, the NRQCD prediction 
on the branching ratio of i/'(3770) — > non - DD which is calculated up to NLO, cannot coincide with the data of BES [2<^. 
At the energy range, the FSI obviously is significant and this allegation has been confirmed by many earlier phenomenological 
studies on other processes. When the FSI effects are taken into account, the discrepancy between theoretical prediction and data 
is significantly alleviated, even though not sufficient. Considering the rather large error range in measurements of both inclusive 
non - DD decay of i^(3770) and the exclusive mode i/'(3770) — > pn, one would still be able to obtain a value for the parameters 
a which does not conflict with the data, by which the theorical prediction and data might be consistent. The more accurate 
mesurements which will be conducted in the future will provide more information which can help to make a definite conclusion 
if the FSI indeed solves the "puzzle" or not. Secondly, our result shows that the FSI can make significant contribution to all the 
channels of i/'(3770) — > pn, K*K, and each of them should be searched in future experiments. Thirdly, no doubt, more accurate 
measurements on i/'(3770) — > non - DD, especially i/'(3770) — > Light Hadrons, are necessary. Thanks to the great improvement 
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of facility and technology of detection at the charm-tau energy region, the BESIII f3^ will provide much more precise data, 
by which we may gain more information. Furthermore, along the other lines more theoretical studies which may involve other 
mechanics, even new physics beyond standard model are badly needed. 
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